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Abstract 



The auxiliary field method, defined through introducing an auxiliary (also called as the 
Hubbard-Stratonovich or the Mean-) field and utilizing a loop-expansion, gives an excel- 
lent result for a wide range of a coupling constant. The analysis is made for Anharmonic- 
Oscillator and Double- Well examples in 0-(a simple integral) and 1- (quantum mechan- 
ics)dimension. It is shown that the result becomes more and more accurate by taking 
a higher loop into account in a weak coupling region, however, it is not the case in a 
strong coupling region. The 2-loop approximation is shown to be still insufficient for the 
Double- Well case in quantum mechanics. 



& ■ 1 Introduction 



In most of actual situations, path integral expressions are given as a non-gaussian form so 
that some approximation is always needed. Apart from perturbative treatment, such as a 
weak(strong) coupling expansion, which can only describe a small (large) coupling region, 
we have sought for other recipes to be able to handle cases for a wider coupling range: the 
variational methodQ has been well known and applied successfully to the polaron problem]^]. 
The method is combined with an optimization technique and has been actively discussed 0. 
Numerical estimation is also possible once expressed in the path integral form: for instance, 
computer simulations produce a lot of fruitful results such as in Lattice QCD[f||, but, in addition 
to the consumption of money (as well as time), there still lacks something to put symmetry 
onto the lattice; chiral symmetry is a well-known example |J. An advantage in path integration, 
contrary to the operator formalism, is that we can easily switch from one variable to other by 
means of some change of variables, which would open a new possibility. The auxiliary field is 
considered as one of these variables, and was introduced into the model by Gross and Neveu||. 

The Gross-Neveu model is a two-dimensional four-fermion model inspired by the work by 
Nambu and Jona-Lasinio[0; 

2 

£ = ^# + y(^) 2 , (1) 
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where ip has an A^-component. After some discussions, they proposed an equivalent Lagrangian, 



O 



(2) 



Here o has no kinetic term and is eliminated by using the equation of motion yielding to the 
original Lagrangian eq.(|l]). In this sense, we call o as an auxiliary field. 

The scenario is more easily understood by means of path integral 0: the partition function 
(in an imaginary temperature) reads 



VipVif) exp 



i I d 2 x I ipiiftip + ^j- (ipijij 



(3) 



Here introducing the auxiliary field o in terms of the Gaussian integrals, such that 



Per exp 



% I d 2 x^ (ex + gipip) 



(4) 



and inserting into eq.(^) we find 

Z = J doVipV^p e ^ d2xC ' = J doV^V^ exp i J d 2 x itfifiip - y - g^pip 



(5) 



Similar techniques are utilized everywhere nowadays also for a boson quartic interaction || 
instead of the four-fermi interaction. The nomenclature for a-field is, therefore, various; the 
mean-field [10], the Hubbard-Stratonovich field [11] in solid state physics. 

In the actual case, we treat the partition function (which can be obtained through 1 1— > it): 



Z(T) = doVipVipexp 



a 



dt J dx I ipiftip + — + gip4>a 



where the anti-periodic boundary condition for the fermi field, if)(T,x) 
understood. We then integrate out the fermion field to find 



(6) 

^(0, x), should be 



2 

Z(T) = J da exp - J d 2 x°— + N In det (0 + go) = J rfae~ 5M 
S[a)s 



(7) 



d 2 x N In det 

2 



g°) 



Since we look for a vacuum with T i— > oo, we should find a constant solution oq in the equation 
of motion, 

5S 



So(x) 







which gives the gap equation; 



o = 2Ng 2 J 



0o 



d 2 k 



(2vr)2 k 2 + (go ) 2 ' 



(8) 



(9) 
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If o"o is non-zero then the dynamically symmetry breaking occurs: this is the end of the usual 
story. Indeed, the recipe is legitimated if the number of fermion species becomes infinite; 
N — > oo. However, it is not the case for most of actual situations: iV is finite or even 1. We 
wish to know "how accurate is it when N = 1?", which is one of the motivation of this work. 

Moreover there is an alternative motivation: performing the WKB approximation in the 
Double- Well potential we must follow instanton calculationsjl2|], which however is very cum- 



bersome as well as tedious. The simpler is the better in any approximation: once introduce an 
auxiliary field we could avoid such troublesome. To clarify these issues, we pick up the case of 
quartic coupling of bosonic field for simplicity. 

The paper is organized as follows: in §11 a simple model calculation is performed for the 
integral expression. Here we realize the importance of the loop-expansion with respect to the 
auxiliary field (variable) and can find a more accurate result is obtained if taking a higher loop 
correction into account when the coupling, g, is small. However, when g goes larger, higher 
loops cannot always improve a situation. We then proceed to the quantum mechanical model 
in §111, where we compare our results with those obtained numerically to find that the 2-loop 
correction gives a 4%-error for 10~ 3 < g 2 < 10 3 except g 2 ~ O(10 _1 ) in the Double-Well case. 
The final section is devoted to a discussion. 



2 Simple (O-dimensional) Model 



The starting point is 



dx 



cxp 



X 

2 



-X 



(10) 




7riun eWVV 



ex 




e(^ 2 )Ii /4 




'ID 



+1 : x > 
-1 : x < 



and l a {x) is the modified Bessel function. There are two cases depending on a sign uj 2 
CASE(i); uj 2 > 0; (O-dimensional Anharmonic-Oscillator), 



W ,.,4//l„2 



2ixg'' 



K 1/4 



(12) 



l + 0(g 2 ) ; 

where K a (x) also is the modified Bessel function. And 
CASE(ii); to 2 < 0; (O-dimensional Double- Well), 

/ 3 ~ y/2 e^ 2 



(13) 
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Here it should be noted that g 2 = is the essential singularity, that is to say, the Double- Well 
case is non-Borel summable. 

Now introduce an auxiliary field such that 



dy 



exp 



-7T [v + w- 



X 



(14) 



so as to erase the x 4 term when inserted into eq . (|10|) , yielding 



dx dy 



2vr V27T 



exp 



1 

r 



r * ( ^ +iW )-v» e -4 = r 

J — oo V Z7T J— i 



dy 



-oo v27T 

We rewrite the final expression as 



exp 



y - - In + 



(15) 



dy 
v/2T 



exp 



S(y) 



a=l 
2 



(16) 



where we have introduced a parameter, a, which must be put unity in the final stage. We call 
a as the loop-expansion parameter. Next, assume the solution of S'(y) = as y Q ; 



S'(y) = y + 



2(uj 2 + igy) 







(17) 



which can be expressed as 



where Q 2 should obey 



Q 2 -u 2 



2Q 2 ' 



Q 2 = uj 2 + igy > ; 



(19) 

since the Gaussian integral of x eq.(^) must exist. Then perform the saddle point method 
around y to give 



dy 
v/2T 



exp 



S(y ) SW(y ) 



a 

S^(y 



Ala 



2a 

<V-Vo)*- 



(y-y T- S{ ^- 



3!a 



-(y-y ) 



(20) 



a=l 



Making a change of variable, y — y i— > y/y/a, we obtain 

o(2) q(3) q (4) 



r _ .-so/a z 00 

' 3 V 2tt 



dyexp 



o(3) 

2 /— °0 3 ^0 4 



e -S /a 

V 2tt 



dy e- 5 o 2) ^/2 h_ ( , 



o(4) 

4! 



o(3) 2 
— y 6 

2(3!) 2 y 



a=l 
2\ 



0(a 2 ) 



(21) 



a=l 



where we have written 



,L-1 



term is called the L-loop term (L = is the tree term.) From eq.([18| 



ft 2 = 

)2 , ,2\ 2 



uj 2 + V^ 4 + ig 2 



S() 



(ft 2 - LO 2 



9 



2g 2 

3 



5, 



(i) 



8ft 4 ' 



5, 



(2) 



1 + 



2ft 2 



O * o(6) 



5!# 6 
2ft 12 



Using these and performing elementary integrals, we obtain 

/ 



exp 



8ft 4 



ft 2 



\ ft 4 + 



1 + 



V 



35g 6 



24(ft 4 + f 



+ 



8 (ft 4 + f 
329c/ 8 



105# 



10 



128 (ft 4 + f) 64 (ft 4 + £ 



+ 0(4-loop) 



Stated as above, we assign 



/ 



(22) 



(23) 



(24) 



(25) 



/tree = exp 



8ft 4 



h- 



loop — tree a 



ft 2 



\ft 4 + f 



'2-loop 



A- 



loop 



3g 4 



8(ft 4 + 



-^3-ioop = exp 



ft 2 



,8ft 4 ; + 

35g 6 



5g e 



24 (ft 4 + i 

Zg A 

8 (ft 4 + f 



+ 



3290 6 



105# 



10 



128 (ft 4 + 



64 (ft 4 + 



(26) 



24 (ft 4 + i 
Let us analyze the individual case: 

• CASE(i); 0-dimensional Anharmonic-Oscillator. Put uj 2 h- > 1 so that eq.(p3|) reads 

„ 2 _ y^+2^+1 



(27) 



We plot the ratio of Jl-1oo P (L=0,1,2,3) to the exact value in Fig.([TJ). (a) shows the case 
of g 2 < 1 and (b) of g 2 > 1. Details are seen in table (a). 
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9 2 


Exact 


Tree 

TVpp / TT!v" 


1-loop 

1 -loon / Ex 


2-loop 
2-1 nnn / Ex 


3-loop 
3-lnrvn / Ex 


1(T 3 


0.9996 


0.9999 
1.00 


0.9996 
1. 


0.9996 
1. 


0.9996 
1. 


10~ 2 


0.9963 


0.9988 
1.00 


0.9963 
1. 


0.9963 
1. 


0.9963 
1. 


10- 1 


0.9685 


0.9881 
1.02 


0.9664 
1 00 

-L .\J\J 


0.9690 
1 00 


0.9683 
1 00 


1 


0.8386 


0.9149 
1.1 


0.8125 
0.97 


0.8541 
1.02 


0.8277 
0.99 


10 


0.5954 


0.7027 
1.18 


0.5484 
0.92 


0.6195 
1.04 


0.5976 
1.00 


10 2 


0.3672 


0.4510 
1.23 


0.3300 
0.90 


0.3817 
1.04 


0.3790 
1.03 


10 3 


0.2131 


0.2656 
1.25 


0.1899 
0.89 


0.2210 
1.037 


0.2222 
1.043 



table(a) 



• CASE (ii); 0- dimensional Double- Well. Put oj 2 ^ — 1 so that eg. (|23|) reads 

= yrw-i m 

2 v ; 

We plot the same ratio as the above in Fig.@. However, in (c) we have omitted the tree 
graph because of the large deviation. Details are shown in table (b). 



9 2 


Exact 


Tree 
Tree/Ex. 


1-loop 
1-loop/Ex. 


2-loop 
2-loop/Ex. 


3-loop 
3-loop/Ex. 


10~ 3 


1.986 xlO 217 


1.035 xl0 2iy 
52.1 


2.313 xl0^ iY 
1.17 


1.930 xl(P Y 
0.97 


1.961 xl0 2iY 
0.98 


io- 2 


7.360 xlO 21 


1.210 xlO 23 
16.44 


8.495 xl0 2i 
1.15 


7.162 xl0 2i 
0.97 


7.270 xl0 2i 
0.98 


HT 1 


2.208 xlO 2 


1.107xl0 3 
5.01 


2.311 xlO 2 
1.05 


2.112xl0 2 
0.96 


2.156 xlO 2 
0.98 


1 


2.350 


4.202 
1.79 


1.932 
0.82 


2.155 
0.92 


2.411 
1.03 


10 


0.8074 


1.103 
1.37 


0.6897 
0.85 


0.8137 
1.01 


0.8768 
1.09 


10 2 


0.4040 


0.5196 
1.29 


0.3542 
0.88 


0.4159 
1.03 


0.4295 
1.06 


10 3 


0.2196 


0.2777 
1.26 


0.1942 
0.88 


0.2270 
1.03 


0.2312 
1.05 



table(b) 



From the figures, ([1]) and ([|), it should be noted that the higher loop corrections improve 
the result all the time when g 2 < 1 but not in the strong coupling region as is seen from the 
graphs (b) and (d). Details for the values are listed in the tables (a) and (b). This fact implies 
that the loop-expansion is merely an asymptotic expansion. In the Anharmonic Oscillator case, 
especially the result is satisfactory: the 2-loop result gives a 4%-error for 1CT 3 < g 2 < 10 3 . In 
the Double- Well case, the 3-loop spoils the result at 1 < g 2 but gives a better result at g 2 < 1. 
However, it is still remarkable that the error, under the 2-loop, remains within ~ 8% for a huge 
coupling region, 10~ 3 < g 2 < 10 3 . 

The essential role of the loop-expansion should finally be remarked: if we stop at the g A 
term in the 2- or 3-loop expression eq. (|26|) , the result deviates far away from the true value. 
Therefore we must abandon the coupling constant expansion in the auxiliary field method. 



3 The Quantum Mechanical Model 



Encouraged by the foregoing results, in this section we analyze the quantum mechanical model: 



H — x H x . 

2 2 8 



(29) 



Here again cases are classified into (i) u 2 > 0; Anharmonic-Oscillator. And (ii) uj 2 < 0; Double- 
Well. The partition function is given by 

Z(T) = Tr e~ TH 
Vxexp 



[ (x 2 u 2 2 g 2 4 

at x H x 

\ 2 2 8 



(30) 



x(T)=x(0) 



where x(T) = x(0) designates the periodic boundary condition. Here and hereafter we put 
Ti i — ► 1 and use the continuous representation, 



N dx- 
Vx = lim TT 



\ - « f = i y/2irAt 



At = 



T 

N 



(31) 



Introducing an auxiliary field in terms of the Gaussian identity, 



J Vyexp ~ J Q dt \[y + 



igx 



2\ 2 ' 



so as to erase the quartic term, we obtain 



dt 



' x 2 



+ [uj + igy 



ar y* 
— + — 
2 2 



Z(T) = JvxVyexp -J 

which becomes, after the integration with respect to x, to 

r r r T v 2 1 ( d 2 N 
Z(T) = JVyexp -J dt y - + - In det I -— +u 2 + igy 



(32) 



(33) 



Vy exp I 



(34) 



a=l 
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where 



%] 



r T v 2 l ( d 2 \ 



(35) 



and again the loop-expansion parameter, a, has been introduced. 

Write a solution of the equation of motion, S'[y] = 5S[y]/5y(t) = 0, y (t), giving the gap 
equation; 

(36) 

(37) 
(38) 



which can be rewritten as 



where 



2/oW + fG(M) = 0, 



Q(t) 2 -u 2 = 9 —G{t,t) 



Q(tf = u 2 + igy (t) 
Here the Green's function, G(t, t'), obeys 

/ d 2 ^ 



dt 2 



Again it should be noted that 



+ VL{tf }G(t,t') = 5(t-t') 



fl(t) 2 > , 



due to the existence of the Gaussian integration of x in eq . (|33|) . 
Now S[y] is expanded around y (or fl(t) 2 ) such that 

%] = S + \{y - y ) 2 • 4 2) + i(y - y ) 3 ■ 4 3) + 



where abbreviations 



S, 



(n) 



S n S 



6y(t 1 )6y(t 2 )---5y(t r 



y=y 



and 



dt x dt 2 ■ ■ ■ dt n (y - y ){ti){y - y ){t 2 ) ■•■{y- y ){t n ) 



5 n S 



y=y 



Sy(ti)6y(t 2 ) ■ ■■Syitr, 

have been adopted. Shifting and scaling the integration variables as before, we obtain 
Z{T) = e' So/a JVyexp 



1 A -1 2 r y 3 c(3) y 4 c( 4 ) 

-A -y -^a — -Sy - a — -Sy 



a=l 



where we have written i— > A 1 which reads explicitly 



6 2 S 



Sy(h)5y(t 2 



y=y 



6{t 1 -t 2 ) + ^G{t u t 2 )G{t2,t l ) . 



(39) 
(40) 



(41) 



(42) 



(43) 



(44) 



(45) 
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Moreover, 

sf = -i^j^i' *2)G(t 2 , * 3 )G(* 3 , *i) + G(t h t 3 )G(t 3 , t 2 )G(t 2 , h) | , 
s£ 4) = -9 4 ^G(h, t 2 )G(t 2 , t 3 )G(t 3 , U)G(U, U) + G(h, t 2 )G(t 2 , U)G(U, t 3 )G(t 3 , h) (46) 

+G(t 1 ,t 3 )G(t 3 ,t 2 )G(t 2l t i )G(U,t 1 )} . 

From these we have 

Z(T) Tree = exp (S ) , 

Z ( T )i-ioo P = ex P (-^o - 2 ln det A ) > (47) 
Z(T) 2 _ loop = exp - g lndet a) x [l + (2-loop Graphs) 

where the 2-loop graphs are formally given by Fig.(|3|). 

Therefore the rest of the work is to fix the form of the Green's function eq. fl39|) and find 
the solution y (t) of the gap equation eq. (|37|) . In this paper we confine ourselves to a time- 
independent solution: those quantities obtained so far turn into the overlined ones; 

y (t) i — ► : constant ; Q(t) 2 i— > Q 2 : constant . (48) 

The Green's function eq.(|30) can be obtained explicitly 



1 oo i2irr(t-t')/T 

G(t,t' ; n) = - £ , 2 , 2 _ 2 



- t') cosh n f — - t + t'^ + 6{t' - t) cosh n(- + t-t' 



2Q sinh ^ 



(49) 



where we have taken the periodic boundary condition into account. 
Now Q is the solution of the gap equation; 



W = - u; 2 = y *; H) = |= coth ^— J . (50) 
When T — > large, Q can be expressed as 

n = n + n ie - naT + n 2e - 2n ° T + ■ ■ ■ , (51) 

where Qq is the solution of the third degree equation; 

njj - uj 2 n = ?L . (52) 
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(In order to calculate the energy of the first excited state, Q± must be known, which is easily 
obtained to be 



g 1 



(53) 



2(3fi§ + u; 2 ) ' 

In this paper, however, only the ground state energy is considered.) Other overlined quantities 
are found straightforwardly; especially 

-l 



A(M') = [ 



PS 



\ Sy(t)6y(t') 



5(t - 1') 



y=y , 



2H 



G(t,t';Cl 



where 



2tt 



The tree part in eq.(^) becomes 



Z ( T ) T rcc = exp(-S, 



with 



■2 i oo 

5o = T^- + i £ In 



2 
T 



2irr 
~T~ 



+ tt 2 



2g 2 



[Q — uj J + In sinh I j + (^-independent part 



(54) 



(55) 



(56) 



(57) 



where the gap equation eq.(|50|) has been utilized to the first term in the final expression. In 
the 1-loop part of eq. fl47f) , we should know lndet A(t, t') which is 



lndet A(t,t') 



1 oo 

5? < ln 



£ M 1 + 



g 2 T 2 



1 



r 2 + 



'my 

27T ) 



-In 



r 2 + 



m 



m/% 

2 



In 



/sinh {VLT/2) 
\ sinh(OT) 



(5f 



As for the 2-loop part, the (non-local) vertices and Sq^ are now expressed as in Fig.@. 
Accordingly, the 2-loop part is written as in Fig. ([5]): 

From the graphs, we should note that there need, in the ordinary sense, the 3- and 4-loop 

calculations in the 2-loop of the auxiliary field, since our vertices, and S^, are non-local. 
Due to this complexity, we confine ourselves to the case that T \— > oo, that is, to the ground 
state. Write the Fourier transformed G and A as 



G{t } t j VLq 
A(t,t') = 



00 dk 



-ik(t-t')_ 



oo 27T 

00 dk 



k 2 + n 2 



00 dk 



ik ^G (k) 



lk(t _ tl) k 2 + 4^ 2 ^ f °o dk 



2vr 



where 



k 2 + Q 2 



Q 2 = AQ 2 + 9 



oo 27T 

(i lfc( '-*' ) A (k) 



(59) 



2vr 



2fi r 



(60) 
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since Q is now Q under T — > oo. With these, each graph, (a) ~ (d), can be expressed and 
calculated elementally as follows: 



j J A (p)A (k)G (l)G (l + fc)G (Z + p)G (l + k + p) 

9 4 P+10 
64fijj p 2 (p + l)(p + 2) ' 

# 4 /■ d/dpdA; 2 # 4 p + 10 



(c) = -^Ao(O) 



(2tt) 2 
g 6 f dl dp dk dq 



A (p)G (l) 2 G (l+p) 



I (p + 6) s 
128Q 8 p 4 (p + 2) 2 



(d) = _^|^~ Ao(p)Ao(?)Ao(p + ?) 



(2tt)< 

xG (l)G (l + p)G (l +p + q)G (k)G (k + q)G (k + q + p) 
g 6 (p 2 + 8p + 4) 



64^p4(p + i)(p + 2 )2 • 
Here we have introduced a parameter, 

P 

The result for the ground state energy, 



w - 



is therefore 



77'tree 


^0 

2 








771 1 — lOOp 
H, Q — 


^0 

2 


(P 


-1)- 


9 2 
32Q 


7-1 2— loop 
'-n 


2 


(P 


-1)- 


9 2 
32Q 



2 ' 


4 



g A (p + 3)(p+10) 



+ 



32ft 2 , 64^ p2( p+1 ) (p + 2 ) 
g 6 p 3 + 15p 2 + 64p + 44 



128Q 8 p 4 (p+l)(p + 2) 2 
Let us analyze the individual case: 



(61) 



0>) - t/w AoW ^^ (62) 



(63) 



(64) 



(65) 



E = -hm ^lnZ(T), (66) 

1 -^oo ± 



(67) 
(68) 



(69) 
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CASE(i); Anharmonic-Oscillator. The solution of eq.fl52"|) is given[[H|] by 



2lu 

7s 



cos 



cos 



-i 



/ 3y/3V 



q 2 to 3 
< — < 



8 V 64 27 \ 



8 " 3^ 



cj 3 o 2 

< y 



64 27 3^ ~ 8 



(70) 



Putting u 2 h 1 we calculate the ratio of E loop (L=0,l,2) to the exact numerical value 
in table (c). 



9 2 /8 


Exact 


Tree 
Tree/ Ex. 


1-loop 
14oop/ Ex. 


24oop 
2-loop / Ex. 


1(T 3 


0.50075 


0.50025 
0.9990 


0.50075 
1. 


0.50075 
1. 


io- 2 


0.50726 


0.50248 
0.9906 


0.50737 
1.0002 


0.50725 
0.9999 


to- 1 


0.55915 


0.52290 
0.9352 


0.56435 
1.009 


0.55775 
0.9975 


1 


0.80377 


0.65268 
0.8038 


0.85522 
1.0640 


0.78548 
0.9772 


10 


1.5050 


1.1080 
0.7362 


1.6729 
1.1116 


1.4738 
0.9793 


10 2 


3.1314 


2.2356 
0.7139 


3.5280 
1.1267 


3.0865 
0.9857 


10 3 


6.6942 


4.7445 
0.7088 


7.5659 
1.1302 


6.6112 
0.9876 



table(c) 

CASE(ii); Double- Well. The solution of eq.(H) is 



\ 



64 + 27 



+ 



\ 



64 + ~W 



(71) 



Putting ui 2 



-1 we again compare the result to the exact numerical value in table (d). 
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9 2 /8 


Exact 


Tree 

Trpp / Th!v" 


1-loop 

1-1 hot) / T^/v 


2-loop 
2-1 nnn / Ex 

_ lKJKJyJ 1 1 J -A . 


icr 3 


-61.794 


-62.500 
1.0114 


-61.794 
1. 


-61.794 
1. 


10~ 2 


-5.5532 


-6.245 
1.1246 


-5.5575 
1.0008 


-5.5541 
1.0002 


10- 1 


-0.15413 


-0.57593 
3 7368 


-0.02326 
1 509 


-0.08479 
5501 


1 


0.51478 


0.25 
0.4856 


0.66421 
1.2903 


0.41605 
8082 


10 


1.3716 


0.92366 
0.6734 


1.5839 
1.1548 


1.2112 
0.8830 


10 2 


3.0695 


2.1501 
0.7005 


3.4867 
1.1359 


2.7543 
0.8973 


10 3 


6.6655 


4.7048 
0.7059 


7.5467 
1.1322 


6.0004 
0.9002 



table(d) 



From the above tables, in case (i) the auxiliary field method can fit the data within a 13%- 
error under the 1-loop and a 3%-error under the 2-loop, which is considered to be excellent. 
For the Double- Well case, the method gives us ~ 10%- error except the region, 0(1CT 2 ) < g 2 < 
0(1), where as was in the 0-dimensional case there might need the 3-loop correction to improve 
the result. Apart from this, it would be still a good approximation for a huge coupling region. 

4 Discussion 

The auxiliary field combined with the loop-expansion can give an excellent result for a huge 
coupling region, O(10~ 3 ) < g 2 < O(10 3 ), even a component of the original variable is single. 
However, in the quantum Double- Well case, there need higher order corrections than the 2-loop 
between O(10~ 2 ) < g 2 < 0(1). A maximum deviation, in the ground state energy of 2-loop, 
reaches 18 times to the exact value with the wrong sign at g 2 ~ 0.15. We have calculated the 
first excited energy E t up to the 1-loop, 

A rr — p 1 - 10 ^ cii-ioop _ 2 ^o ( i 9 2 3g 2 j 1 A 

AE = El ~ E ° -MT^V'M ^miMT^ 2 ) (72) 

and found a level crossing around these regions. Apparently the approximation is broken down 
there. However, it is cumbersome to go beyond the 1-loop in quantum field theory as well 
as quantum mechanics. The approximation scheme should be simple and transparent. We 
therefore look for another solution rather than a time-independent solution, that is, we must 
solve eq. (|39|) more carefully. Indeed, the structure of the dominant contribution to path 
integral has recent been clarified by means of such as the valley method [p^]. With these in 
mind the work is in progress. 

As for applications, the recipe is applicable almost to any situation. Our interest is the 
dynamical structure of QCD that is recently revealed in terms of a profound consideration into 
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gauge invariance by Lavelle and McMullan et al.[15|, [16] , for example. It is thus tempting to 
introduce this method into QCD, which is also our future program. 
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3-1 oop 



— Q 
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l-loop 



(a) 



(b) 



Figure 1: 0- dimensional Anharmonic Case: (a); g 2 < 1. (b); g 2 > 1. Dotted line: /Tree/Exact. 
Dashed line: ii_i oop /Exact. Solid line: /2-ioop /Exact. Dash-dotted line: /3_i oop / Exact. 




(c) (d) 

Figure 2: 0-dimensional Double-Well Case: (c); g 2 < 1. (d); g 2 > 1. Dotted line: /Tree/Exact: 

this is omitted in (c); because of a large deviation. Dashed line: Ii-\ oop / Exact. Solid line: 
h -loop /Exact. Dash-dotted line: /3_i o P /Exact. 
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Figure 3: Formal 2-loop Graphs: A denotes the propagator of the auxiliary field. 




Figure 4: Vertices for the constant classical solution: the solid line represents G\t, t';fl) and 
the dot represents the vertex which should be attached by the double-lined A(t,t'). 
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(a) (b) ( C ) (d) 



Figure 5: The 2-loop Graphs for the constant classical solution: the upper numbers in the 
figures represent those of multiplicity. 
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